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RELATIVISTIC  FOCUSING  AND  BEAT  WAVE  PHASE  VELOCITY 
CONTROL  IN  THE  PLASMA  BEAT  WAVE  ACCELERATOR 

Introduction 

Recently  there  hsis  been  much  interest  in  pl2tsma  based  accelerator  schemes,  such  eis 
the  plasma  beat  wave  accelerator  (PBWA),^~’  for  producing  ultra-high  energy  electrons. 
This  has  led  to  a  renewed  interest  in  the  study  of  the  propagation  of  intense  radiation  beams 
through  a  plasma.^”^®  In  the  PBWA  two  colinear  radiation  beams  of  frequencies  u)i,u>2 
are  incident  on  a  umform  plasma.  By  appropriately  choosing  the  difference  in  the  lasci 
frequencies  to  be  equal  to  the  electron  plasma  frequency  u>p,  Auf  =  u>i  —  u>2  =  u/p,  where 
uip/uij  <<  1,  it  is  possible  for  the  radiation  beat  wave  to  resonantly  drive  large  amplitude 
electron  plasma  waves.  In  the  ideal  wave  breaking  limit, the  maximum  accelerating 
electric  field  Em  is  given  by  Em  =  (meC*/e)wp/c  ^  .97^/np  eV/cm  where  Up  is  the  plasma 
density  in  cm“®.  For  example,  Up  =  1.6  x  10**  cm~*  gives  Em  —  120  MeV/cm  which 
implies  that  an  electron  could  be  accelerated  to  1.2  TeV  in  100  meters. 

To  realize  such  an  acceleration  scheme  it  is  necessary  that  t)  the  radiation  beams 
propagate  at  high  intensity  over  distances  large  compsired  to  the  Rayleigh  length  zr  = 
u;rj/2c,  where  r,  is  the  radiation  spot  size,  and  that  it)  phase  resonance  between  the 
accelerating  electrons  and  the  plasma  wave  be  maintained  over  an  equally  large  distance. 
In  vacuum,  radiation  diffracts  over  distances  on  the  order  of  cr,  which  can  be  relatively 
short.  Hence,  in  order  to  maintain  high  intensity  beams  it  is  necessary  to  rely  on  focusing 
enhancement  from  the  plasma.  In  the  PBWA  the  phase  velocity  of  the  plasma  wave  is 
equal  to  the  phase  velocity  of  the  radiation  beat  wave  which,  in  the  1-D  limit,  is  given  by 
Vp/c  =  A(jj/Ak  ~  1  —  ci;p/2wf,  where  Ak  =  ki  —  k2  is  the  difference  in  the  wave  numbers 
of  the  two  beams.  Since  the  velocity  of  an  ultra-relativistic  electron  is  approximately  the 
speed  of  light,  the  electrons  out  run  the  plasma  wave  and  become  “detuned”  in  a  length*® 
Ld  —  where  Xp  =  2vc/u)p.  For  wi/u;p  =  25  and  Up  =  1.6  x  10*®  cm~*,  this  gives 

Ld  —  16  cm  and  a  maximum  electron  energy  gain  of  A€  ~  EmLd  —  2  GeV.  In  order  to 
increase  the  energy  gain  beyond  this  detuning  limit,  it  is  necessary  to  increase  the  pliase 
velocity  of  the  plasma  beat  wave. 

Tliis  paper  addresses  the  two  points  mentioned  above  concerning  the  realization  of 
the  PBWA.  As  is  shown  below,  matched  beam  solutions  are  possible  in  which  the  two 
radiation  beams  propagate  with  constant  spot  sizes  provided  the  radiation  is  of  sufficiently 
high  power.  This  adlows  the  radiation  beams  to  propagate  over  distances  larger  than 
the  Rayleigh  length  while  maintaining  their  high  intensities.  For  example,  two  radiation 
beams  with  equal  spot  sizes  are  matched  provided  the  power  in  each  beam  is  P  =  P,.  3. 

Manuscript  approved  July  7.  1988. 
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where  Pc  ~  17  x  10*u;^/u;p  W  is  the  power  threshold  foi  relativistic  focusing  of  a  single 
radiation  beam  in  a  plasma.^ In  addition,  by  including  finite  radial  beam  profiles  along 
with  relativistic  focusing,  the  phase  velocity  of  the  beat  wave  can  be  tuned  to  the  speed 
of  light.  This  is  accomplished  by  appropriately  choosing  the  initial  spot  sizes  and  powers 
of  the  radiation  beams.  Hence,  phase  resonance  between  the  electron  and  beat  wave  can 
be  maintained  beyond  the  1-D  detuning  length  Lj,  and,  consequently,  substantially  higher 
electron  energies  can  be  achieved.  Figure  1  shows  schematically  the  propagation  of  two 
matched  radiation  beams  through  a  plasma  with  the  resulting  beat  wave  phase  velocity 
equal  to  the  speed  of  light. 

Focusing  of  radiation  beams  in  a  plasma  occurs  through  the  combined  effects  of  rela¬ 
tivistic,  ponderomotive  and  thermal  self-focusing.^”^®  Typically  these  processes  occur  on 
widely  separate  time  scales.  Relativistic  focusing*”*  occurs  on  the  shortest  time  scale, 
rR  1/w,  which  is  the  time  scale  at  which  the  electrons  respond  to  the  radiation  field. 
Ponderomotive  focusing*”^^  depends  on  the  expulsion  of  ions  from  the  radiation  chan¬ 
nel  and  thus  occurs  on  a  time  scale  given  roughly  by  rp  ~  where  C,  is  the  ion 

acoustic  speed.  Thermal  focusing^®’^®  relies  on  heating  of  the  plasma  by  the  radiation 
beam  and  typically  occurs  on  an  even  longer  scale.  This  paper  is  concerned  with  rela¬ 
tivistic  focusing,  and  hence  the  analysis  is  applicable  to  lasers  with  pulse  lengtlis  in 
the  range  tr  <  ri  <  rp,  which  is  the  region  of  interest  for  the  PBWA.  Physically,  rela¬ 
tivistic  focusing  eurises  solely  from  the  relativistic  electron  quiver  velocity,  r,  =  caj./7j., 
in  the  combined  radiation  field.  Here  aj.  =  eAi^lmc'  is  the  normalized  radiation  vector 
potential  and  7x  =  relativistic  gamma  factor  for  an  electron  in  a  helically 

polarized  radiation  field.  The  focusing  mechanism  for  a  single  beam  is  that  a  radiation 
profile  peaked  on  axis  leads  to  an  index  of  refraction  profile,  n  ~  1  —  (u;j,/u>)®/27x,  which 
has  a  minimum  on  axis.  The  radiation  beam,  therefore,  focuses  along  the  axis.’^  When 
the  radiation  power  is  greater  than  the  critical  power  Pc  for  relativistic  self-focusing,  it  is 
possible  for  the  envelope  of  a  single  radiation  beam  to  propagate  at  a  constant  spot  size.' 
For  two  colinear  beams,  however,  the  situation  is  more  complicated  due  to  the  coupling 
of  one  beam  to  the  other  through  the  relativistic  gamma  factor.  The  analysis  presented 
below  indicates  that  matched  beam  propagation  for  two  beams  is  only  possil)le  for  a  finite 
range  of  the  parameter  iZ,  such  that  l/(\/8  —  !)</?<  v'8  -  1.  where  R  =  (c,i/r,2)®  is 
the  square  of  the  ratio  of  the  spot  sizes  of  the  two  beams.  The  radiation  power  required 
to  obtain  matched  beam  propagation  for  two  beams  is  near  that  for  a  single  beam.  Pr. 

Control  of  the  beat  wave  phase  velocity  is  most  easily  understood  by  the  followinii 
heuristic  argument.  The  finite  radial  extent  of  the  radiation  beams  gives  rise  to  a  small 
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effective  perpendicular  wave  number  kx-  Here  kx  is  a  function  not  only  of  the  spot  size 
but  also  of  the  power  due  to  the  relativistic  focusing  effects.  The  existence  of  kx  gives 
rise  to  an  effective  parallel  wave  number  given  by  fc||  ~  (1  —  u)p/2a;^  -  c~k\l2ijj‘^)(jj j c. 
Hence,  the  parallel  phase  velocity  of  the  beat  wave  is  now  given  by  Vp/c  ~  1  -  u;^/2u;^  + 
—  k\j2)c^ /2uu)p.  By  appropriately  choosing  the  initial  spot  sizes  and  powers  of  the 
two  radiation  beams,  it  is  possible  to  have  the  last  term  in  the  expression  for  Vp/c  cancel 
the  second  term  thus  providing  Vp  =  c. 

Analysis  of  Radiation  Focusing  and  Beat  Wave  Phase  Velocity  Control 

The  analysis  starts  with  the  wave  equation  for  the  vector  potentiad  of  the  ccmhi!i'*d 
radiation  field,  (V*  —  ldi‘)Ax  —  — (4ir/c)Jj.,  where  Jx  is  the  transverse  current 

density.  In  order  to  study  the  effects  of  relativistic  focusing  alone,  only  the  current  resulting 
from  the  electron  quiver  motion  is  needed,  Jx  —  —eripVq.  The  wave  equation  is  then  given 
by  ^ 

(v^  -  (l  +  |ai|^  +  |a2p -h2|ail|a2|cosA$)  ,  (1) 

where  ax  =  nj  +02"  Throughout  the  following,  a  subscript  1  refers  to  the  radiation  beam 
of  frequency  Wi,  and  a  subscript  2  refers  to  the  radiation  beam  of  frequency  W2-  The  factor 
within  the  square  root  is  the  relativistic  gamma  factor  7x>  assuming  helically  polarized 
radiation.  Here,  A#  =  # j  — $2  *s  phase  of  the  beat  wave  and  $1,2  =  ^1,2-  -‘*’1,2 ^  -'■<^1.21 
where  <^1,2  is  the  slowly  evolving  phase  of  the  radiation  field. 

In  order  to  examine  the  qu2ditative  diffractive  properties  of  the  radiation  beams,  it  is 
helpful  to  consider  the  index  of  refraction  ni,2  of  each  beam.  The  approximate  index  of 
refraction  associated  with  each  beam  is  obtained  in  the  following  manner.  First,  the  1-D 
limit  of  the  left  hand  side  of  Eq.  (1)  is  taken,  assuming  a  ~  exp(t^).  Next,  Eq.  (1)  is 
divided  by  the  phase  factor  exp(z$)  of  either  beam  1  or  2  and  then  averaged  over  a  period 
of  the  beat  phase  A$.  In  the  mildly  relativistic  limit  |ai,2p  <  1,  the  index  of  refraction 
for  each  beam  is  given  by 


ni  =  A:ic/u;i  =  1  -  (w'/2u/{)(l  -  lail*/2  -  |o2|‘),  (2a) 

712  =  k2c/(jj2  =  1  —  (u?p/2u;2)(l  —  |ai|‘  —  |a2i"/2).  (26) 

In  the  above  expressions,  the  first  term  (the  unity)  represents  vacuum  diffraction  wliile  the 
second  term  is  the  contribution  from  the  ambient  plasma.  The  remaining  terms  represent 
focusing  from  the  radiation  fields.  More  specifically,  a  term  proportional  to  ( |ai  ('  +  lao  1' )  2 
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results  from  the  individual  contributions  of  beam  1  and  2  to  the  relativistic  factor  7x, 
while  the  remaining  term  proportional  to  |ai,2p/2  results  from  the  contribution  of  the 
beat  wave  to  7x.  Radiation  focusing  occurs  when  dn/dr  <  0.  Hence,  the  contribittion 
from  the  radiation  terms  to  ni,2  provide  focusing  for  radiation  profiles  peaked  on  axis.  For 
sufficiently  high  power,  these  focusing  terms  dominate  the  vacuum  diffraction  and  provide 
overall  focusing  of  the  radiation  besuns. 

Envelope  equations  describing  the  evolution  of  the  spot  size  r,{z)  of  each  beam  are 
derived  by  applying  the  “source-dependent  expansion”  (SDE)^*'^^  to  Eq.  (1).  This  is 
accomplished  by  expanding  the  normalized  vector  potential  ai,2  for  each  beam  into  a 
series  of  Gaussian-Laguerre  polynomials  and  using  orthogonality  properties  to  determine 
their  coefficients.  The  SDE  differs  from  the  typical  vacuum  modal  expansion  in  that  the 
parameters  characterizing  the  Gaussian-Laguerre  polynomials,  such  as  the  width  of  the 
Gaussian,  are  functions  of  z  which  depend  on  the  “source”,  i.e.  the  right  hand  side  of 
Eq.  (1).  Assuming  that  each  beam  is  adequately  described  by  the  lowest  order  Gaussian 
mode  a  =  |aoo|exp[i/?  —  (1  ~  then  the  parameters  |aoo|»  /d,  a  and  t,  are  given 

by  |aoo(2)|  =  <*(2)  =  (a;/4c)drj/dz  along  with  the  following  equations: 


dz2  **  u;?r 


!£Lfi_vr  + 


u;|rj2  [  W,(l  +  R)2)\^ 


d  ^  2c 

— - - 

dz  wirjj 

d  ^  2c 

—  02  - - — 

dz 


,  t2,.2 


n  aw2R{ 


R{l  +  2R) 

(1  +  R0 


1  +  -  W2 

4c- 


AW,(2  +  R)' 
W2{1  +  R)2 


where  the  mildly  relativistic  limit  wais  taken,  |ai,2p  <  1.  Physically,  W  =  (u;paor,o/4c)^  = 
P/ Pc  where  P  is  the  power  in  one  of  the  beams  and  Pc  is  the  critical  power  necessary 
for  relativistic  focusing  of  a  single  beam.'  Here  a®  and  r,o  are  the  initial  amplitude  of 
the  vector  potential  on  axis  and  the  initial  spot  size  of  each  beam.  The  parameter  a  is 
related  to  the  curvature  of  the  radiation  wavefront  and  the  parameter  d  is  important  in 
that  it  represents  a  correction  to  the  parallel  wave  number  on  axis  A;||  =  ^/c  +  dS  jdz.  This 
relation  is  used  below  to  determine  the  beat  wave  phase  velocity  on  axis. 

Equations  (3a)  and  (3b)  describe  the  envelope  evolution  for  each  beam  as  it  propagates 
through  the  plasma.  Setting  the  right-hand  sides  of  Eqs.  (3a)  and  (3b)  equal  to  zero  gives 
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matched  beam  solutions  for  wliich  the  beams  propagate  without  diffracting.  Matched  beam 
solutions  are  obtained  for  values  of  R  in  the  range  l/(\/8  —  1)  <  i?  <  —  1  provided  the 

normalized  power  W  of  each  beiun  is  given  by 

=  [8R\1  +  RY  -  (1  +  R)*][64R^  -  (1  +  R)*]-\  (5a) 

W2  =  [8(1  +  -  (1  +  R)*][64R^  -  (1  +  /?)*]-'.  {5b) 

This  is  illustrated  in  the  following  limits:  For  i?  =  1,  then  =  W2  =  1/3.  As  R  ^ 
l/(>/8  —  1),  then  Wi  —*■  0  amd  W2  — >  1.  As  i?  — ►  v/S  —  1,  then  VTi  — ►  1  and  W2  —  0. 
Hence,  it  is  possible  for  a  beam  close  to  the  critical  power  to  confine  a  second  beam  which 
has  a  smaller  spot  size  and  a  smaller  power. 

.Matched  beam  propagation  occurs  when  R,  VFj  and  W2  are  specified  as  indicated 
above.  For  example,  once  R  is  chosen  in  the  range  1/(n/8  —  l)<i2<v/8-l,  then 
and  W2  are  given  by  Eqs.  (5a)  and  (5b).  The  actual  magnitudes  of  the  spot  sizes  r,i  and 
r,2  are  undetermined  and  only  their  ratio  has  been  specified.  Specifying  a  value  for 
gives  a  value  for  the  radiation  beat  wave  phase  velocity  on  axis  according  to  the  relation 
c/vp  =  1  4-  cAl3'/A(jj,  where  A0'  =  d^i/dz  —  dl32ldz.  Alternatively,  requiring  Vp  —  c 
for  a  given  set  of  matched  beam  parameters  R,  Wi  and  IV2  specifies  r,i.  For  example, 
as  ii  -♦  \/8  -  1,  requiring  Vp  =  c  gives  ~  5u>i/Au>,  where  kp  =  ujp/c.  For  i?  =  1, 

requiring  Vp  =  c  gives  ~  2.  As  i?  —  l/(v^-  1),  it  is  not  possible  to  have  Vp  =  c.  For 
applications  in  the  PBWA,  it  may  be  desirable  to  have  >>  1.  This  implies  that  it 
may  be  desirable  to  choose  a  matched  beam  case  with  i?  >  1.  For  example,  R  =  1.5  gives 
Wi  ~  0.7,  W2  ~  0.1  and  ~  2.6u}i/ Au). 

As  a  final  illustration,  the  above  results  are  applied  to  parameters  similar  to  those  in 
the  UCLA  beat  wave  excitation  experiment,**  where  u/j  ~  2.0  x  10**  sec“*  and  Au;/u;i  tr 
9.7  X  10~^  (wliich  implies  Up  ~  10*^  cm~*).  A  test  electron  with  intial  energy  given  by 
7o  =  50  is  accelerated  by  plasma  waves  generated  in  the  following  two  special  cases:  i)  A 
matched  beam  case  with  the  beat  wave  phase  velocity  tuned  to  the  speed  of  light.  Vp  —  c. 
where  R  =  1.5,  r,i  =  8.3  x  10"*  cm.  Pi  =  1.3  x  10*-  W  and  P.  =  1.6  x  10**  W;  and 
»i)  the  same  parameters  as  case  i)  only  now  the  beat  wave  phase  velocity  is  given  liy  the 
1-D  limit,  Vp/c  ~  1  —  u7p/2u;j,  and  the  radiation  beams  are  assumed  to  undergo  vacuum 
Rayleigh  diffraction,  r,  =  r,o(l  +  The  results  of  caise  ii)  are  shown  in  Fig.  2 

and  the  results  of  case  t)  are  shown  in  Fig.  3.  Figure  2  indicates  that  the  test  electron 
outruns  the  plasma  wave  and  begins  to  be  deaccelerated  after  approximately  0.5  cm  with 
a  maximum  energy  gain  of  A7  ~  210.  In  Fig.  3,  however,  phase  resonance  between  the 
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electron  and  beat  wave  is  maintained  which  allows  an  energy  gain  of  Aj  ~  6500  in  H  cm. 
This  energy  gzun  continues  in  a  linear  fashion  until  it  becomes  limited  by  some  non-ideal 
effect  such  as  pump  depletion.** 

Discussion 

In  summary,  it  has  been  shown  that  two  colinear,  short-pulse,  Gaussian  radiation 
beams  can  propagate  through  a  uniform  plasma  without  diffracting  due  to  relativistic 
focusing.  This  occurs  for  values  of  Jl  in  the  range  l/(\/8  —  1)  <  <  \/8  —  1,  provided  Wi 

and  W2  are  specified  according  to  Eqs.  (5a)  and  (5b).  In  addition,  it  is  possible  to  tune 
the  phase  velocity  of  the  radiation  beat  wave  to  the  speed  of  light  for  cases  where  iZ  >  1. 
This  is  accomplished  by  appropriately  choosing  r^i.  In  an  actual  PBWA,  Au)  =  u?p  and  the 
envelope  behavior  of  the  radiation  beams  becomes  more  complicated  due  to  the  presence 
of  large  amplitude  resonantly  driven  plasma  waves. However,  the  analysis  presented  here 
remains  valid  for  the  front  of  the  radiation  pulse  (the  first  several  plasma  wavelengths) 
where  the  amplitude  of  the  plasma  wave  remains  small.  In  the  small  amplitude  limit,  the 
phase  velocity  of  the  plasma  wave  is  equal  to  that  of  the  radiation  beat  wave.^*  .4ssuming 
that  the  phase  velocity  of  the  plasma  wave  remains  fixed  to  its  initial  value,  then  the  above 
analysis  indicates  that  it  is  possible  to  tune  this  phase  velocity  to  the  speed  of  light.  This 
implies  that  phase  detuning  between  the  plasma  wave  and  the  electrons  can.  in  principle, 
he  avoided  which  results  in  »  substantially  higher  energy  (zain  in  the  PBWA. 
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Energy  ( MeV ) 


FIG.  2.  Test  electron  acceleration  with  the  same  parameters  as  Fig.  3  except  the  phase 
velocity  is  given  by  the  1-D  limit  and  the  radiation  beams  undergo  vacuum  Rayleigh 
diffraction. 
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Energy  ( GeV ) 


FIG.  3.  Test  electron  acceleration  for  a  matched  beam  case  with  Vp  =  c,  where  R  =  1.'. 
r,i  =  8.3  X  10“*  cm,  Pi  =  1.3  x  10^^  W  and  P2  =  1.6  x  10^^  W. 
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